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ABSTRACT 



Context. The stellar mass function is a probe for a potential dependence of star formation on the environment. Only a few young 
clusters are known to reside within the Central Molecular Zone and can serve as testbeds for star formation under the extreme 
conditions in this region. 

Aims. We determine the present-day mass function of the Quintuplet cluster, a young massive cluster in the vicinity of the Galactic 
centre. 

Methods. We use two epochs of high resolution near infrared imaging data obtained with NAOS/CONICA at the ESO VLT to 
measure the individual proper motions of stars in the Quintuplet cluster in the cluster reference frame. An unbiased sample of cluster 
members within a radius of 0.5 pc from the cluster centre was established based on their common motion with respect to the field and 
a subsequent colour-cut. Initial stellar masses were inferred from four isochrones covering ages from 3 to 5 Myr and two sets of stellar 
evolution models. For each isochrone the present-day mass function of stars was determined for the full sample of main sequence 
cluster members using an equal number binning scheme. 

Results. We find the slope of the present-day mass function in the central part of the Quintuplet cluster to be ff = -1.66 ± 0.14 for 
an approximate mass range from 5 to 40 Mq, which is significantly flatter than the Salpeter slope of a = -2.35. The flattening of 
the present-day mass function may be caused by rapid dynamical evolution of the cluster in the strong Galactic centre tidal field. The 
derived mass function slope is compared to the values found in other young massive clusters in the Galaxy. 

Key words. Galaxy: center - Galaxy: open clusters and associations: individual: Quintuplet cluster - stars: luminosity function, mass 
function - stars: early-type - infrared: stars - instrumentation: adaptive optics 



1. Introduction 

The Quintuplet cluster is one of only three young, massive clus- 
ters known within the central molecular zone (CMZ) with pro- 
jected distances of less than 30 pc to Sagittarius A*, the super- 
massive black hole (SMBH) at the centre of the Milky Way. The 
other two clusters are the Arches cluster at a similar location 
as the Quintuplet cluster, and the nuclear star cluster surround- 
ing Sgr A*. These clusters are unique laboratories to study the 
formation and evolution of stars and their host clusters in the 
Galactic Centre (GC) environment. 

The conditions for star formation in the CMZ and the GC 
region are rather extreme in terms of high gas densities, en- 
hanced temperatures, tidal forces exerted by the gravitational 
potential in the inner Galaxy, and strong magnetic fields. These 
conditions were suggested to favour the formation of high mass 
stars as compared to the more moderat e spiral arm environments 
(iMorrisI I 19931 1m orris & Serabvnl[T996i) . An overpopulation of 
high mass stars may be evidenced in a flattened initial mass 
function (IMF) in GC star clusters. The young massive clus- 
ters are ideal candidates to search for such a deviation from the 
Galactic field IMF. Their youth ensures that a large fraction of 
the initial population is still present in or near the cluster and 
their high total masses provide coverage of the entire mass func- 



* Based on observations collected at the ESO/VLT under Program ID 
71.C-0344(A) (PI: F. Eisenhauer, retrieved from the ESO archive) and 
Program ID 081.D-0572(B) (PI: W. Brandner). 



tion (MF) up to the highest- mass stars known, such as the Pistol 
star in the Quintuplet cluster (Fig er et al.ll995l:[Figer et al.ll998l: 
lYungelson et al. 2008) . Due to the large number of high mass 
stars these cluster are also well-suited to assess stellar evolution 
scenarios for the most massive stars. 

A direct comparison of the observed present-day mass func- 
tion (PDMF) of the GC young massive clusters with the Galactic 
field IMF is aggravated due to their rapid dynamical evolu- 
tion and dissolution in the GC tidal field. N-body simulations 
of compact massive clusters with m asses < 2 xJ_O^M0 and 
distances from the GC < lOOpc bv iKim et all (l2000h yielded 
dissolution times of less t han 10 Myr. A similar study by 
Portegies Zwart et al. (2002) derived somewhat longer dissolu- 
tion times of up to 55 Myr for a GC distance of 150 pc, but found 
that the spatial density of a young massive cluster drops quickly 
below the background density within a few Myr, rendering older 
clusters indetectable. 

In spite of this difficulty, measuring the PDMF is essential to 
deduce the IMF and to compare star formation in the GC with 
the outcome in less extreme environments. As the cluster centre 
is least affected by tidal stripping, as opposed to the outer clus- 
ter areas, the slope of the PDMF there can be compared to the 
PDMF slopes of spiral arm young massive clusters. In addition, 
tidally stripped stars might still be located close to the cluster 
at these young ages. If these stars can be identified as former 
cluster members, e.g. as they are co-moving with the cluster at 
comparable velocities, the common PDMF of cluster members 
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and former members might still be a reasonable representation of 
the IMF. Most notably the measured PDMF of a GC cluster is an 
indispensable ingredient to constrain dynamical simulations set 
out to reconstruct its dynamical history, its IMF and to recover 
possible IMF variations. 

The name of the Quintuplet cluster arises from a "quintuplet" 
of bright near-inf rared (NIR) sources (01 -04,09), which was 
first observed by iNapata et alj (Il990h and lOkuda et al.1 d 19901) 
by resolvi ng one or two merged sources present in earher sur- 
veys (e.g. iKobavashi et an il983 ). Shortly afterwards, the num- 
ber of resolved clu ster stars was increased to 15 star s (01-015, 
iGlass et aL I dun, the "0"-label was introduced in iFiger et al] 
(Il999bl) '). Until today, 21 Wolf-Rayet (WR) stars (14 WC, 
7 WN stars, see [va n der Hucht 20 06i and referen ces therein; 
iMauerhan etaTIIIoTobt iLiermann etaLll2009l 1201 Ol) . two lumi- 
nous blue variables (LBVs) and 93 OB stars were spectro - 
scopically identified dFiger et al.lll999bt ILiermann et al. I I2OO9I) . 
Hence, the Quintuplet hosts almost a quart er of the 92 WR stars 
in the Galactic centre region (Mauerhan e t alJl2010ah . The age 
of the cluster was derived b ased on spe ctral types of likely clus- 
ter stars to be 4 + 1 Myr ( IFiger et al.llT999b.) . The Quintuplet 
cluster is the refore slightly older than the Arches cluster (age: 
2-2.5 Myr, Nai arro et al.l 12004*). with whom it shares a simi- 
lar projected GC distance (30 pc and 26 pc), a similar location 
in the sky and comparable masses on the order of lO'^Mo. As 
the density of the Quintuplet cluster is more than two ord ers of 
magn itude smaller than the density of the Arches cluster dFigeii 
12008"), and it exposes a much more dispersed appearance, it is 
suggestive to regard these two clusters as snapshots at different 
timesteps of massive clusters rapidly dissolving in the GC tidal 
field. 

Due to the strong interstellar extinction for lines of sight to- 
wards the GC and the large distance to these clusters, NIR ob- 
servations with high spatial resolution provided either by ground 
based telescopes with AO correction or space telescopes are re- 
quired to study their stellar populations. As the Quintuplet clus- 
ter is far less compact than the Arches cluster, it is mandatory 
even in the cluster core to apply an effective method to discern 
the cluster stars from field stars in order to determine cluster 
properties in an unbiased way. In this paper, we analyse adap- 
tive optics (AO) observations of the Quintuplet cluster obtained 
with the NAOS-CONICA (NACO) instrument on the ESO Very 
Large Telescope (VLT). A clean sample of cluster stars was es- 
tablished based on the measurement of individual proper mo- 
tions and a subsequent colour-cut. The PDMF of the central part 
of the Quintuplet cluster is derived from this cluster sample for 
the first time. 

In Sect. |2] we present the data and the data reduction. The 
photometric calibration and the determination of the astrometric 
and photometric uncertainties are detailed in Sect. [3] The com- 
pleteness of detected stars is derived from artificial star experi- 
ments in Sect. |4] The proper motion of individual stars is deter- 
mined in Sect. |5j and the proper motion membership sample is 
established. The colour-magnitude diagrams (CMDs) of proper 
motion members and non-members are discussed in Sect. |6l 
and the final cluster sample is derived applying a colour-cut 
and using spect r al iden tifications from the spectral catalogue by 
ILiermann et al.l d2009l) . Initial stellar masses are determined in 
Sect.|7]from four isochrones of different ages and different stellar 
models. The arising PDMFs of the four isochrones are discussed 
in Sect. [8] In Sect. |9] we compare our results to reported mass 
function slopes of other young massive clusters in the Milky 
Way. This paper concludes with a short summary in Sect. [10] 



2. Observations and data reduction 

Two epochs of high-resolution observations of the Quintuplet 
Cluster (a = 17H6'"15\ 6 = -28°49'41", J2000) were ob- 
tained in service mode on July 22-23th 2003 and July 24th 
2008 with the Very Large Telescope (VLT), yielding a time 
baseline of 5.0 yr. In order to achieve the high astrometric ac- 
curacy crucial for this study, the near-infrared imaging camera 
CONICA was used in combination with th e NAOS instrument 
prov iding adaptive optics (AO) correction dLenzen et al J 120031 : 
Rou sset et alJl2003l) . The first epoch was observed with a Santa 
Barbara Research Center (SBRC) InSB Aladdin2 array as detec- 
tor, which was replace d by an InSB Aladdin3 array in May 2004 
dAgeorges et al.l [20071) . All data were taken with the medium- 
resolution camera (S27) with a pixel scale of 0.0 27 T'pixeL^ 
and a field of view of 27.8". The Quintuplet star Q2 dGlass et al.l 
Il990 !) with Ks ~ 6.6 mag served as the natural guide star for 
the infrared wavefront sensor. The data of the first epoch was re- 
trieved from the ESO archive (PI: F. Eisenhauer) and consisted 
of three datasets: two datasets in the H- and /Tj-band with a short 
detector integration time (DIT) of 2.0 s and NDIT (number of in- 
tegrated DITs) of 30 and a further /Tj-band dataset with a longer 
DIT of 20.0 s and NDIT of 2 to expand the range of observed 
magnitudes towards fainter stars. Each of these three datasets 
consists of 16 dithered frames resulting in an effective field of 
view of about 40" x 40". The total integration time at the cen- 
tral part of the observed field with maximum overlap was 16 min 
for the two datasets with the shorter DIT and 1 1 min 40 s for the 
dataset with the longer DIT. The /Tj-band observations of the 
second epoch in 2008 cover the same field of view as the first 
epoch, but were carefully designed to provide high astrometric 
accuracy. The DIT of this dataset was 2.0 s with NDIT of 15, re- 
sulting in a total integration time of 16 min 30 s for the 33 frames 
used for the combined image out of a total of 44 frames. The 
properties of all four datasets are summarized in Table [T] and a 
composite colour image is shown in Fig.[T] 

Standard data reduction was performed using a custom made 
reduction pipeline written in Python, calling a series of self- 
written IDL routines as well as PyRAF tasks. Both sky and sci- 
ence frames were reduced by subtracting the appropriate darks 
and dividing by the twilight flat fields. One master sky was cre- 
ated per dataset with the PyRAF/IRAF task imcombine by deter- 
mining the average of the 2nd to 5th faintest pixel. All brighter 
pixels were rejected to avoid contamination with stellar light in 
the master sky image. This is particularly crucial along the line 
of sight towards the GC, where the high stellar density prohibits 
the choice of star- free sky fields. For the 2008 dataset, the mas- 
ter sky was derived from both the sky and science frames as oth- 
erwise strong stellar residua remained resulting from too small 
dithers of the sky frames. Before the master sky image was sub- 
tracted from a science frame it was scaled to the background 
level of this frame. Hot and dead pixels were detected from out- 
liers during the combination of the dark images and twilight flats 
and written to the master bad pixel mask. For each science frame, 
pixels affected by cosmic ray hits were determined by the IRAF 
task cosmicrays and added to the master bad pixel mask to create 
an individual mask for each frame. 

All data from 2008 except for the twilight flat fields were af- 
fected by the sporadic so-called 50 Hz noise from the Aladdin3 
detector causing horizo ntal stripes, the in tensity and position of 
which vary with time (Lundin et al. "2007, chap. 5.1). To correct 
for this noise, an appropriate routine from the ESO ECLIPSE 
pipeline dDev illard 2001) for the ISAAC instrument was rewrit- 
ten in IDL and adapted to the NACO data. It encompasses the 
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following steps: 1.) Determine the median of each row after re- 
jecting the 40 darkest and the 420 brightest pixels in order to ex- 
clude bad pixels and flux from stars, and store the median values 
in an array. 2.) Smooth the array of median values with a median 
filter with a half-width of 40 pixel. 3.) Subtract the difference 
between the original and the smoothed aiTay of median flux val- 
ues from each column of the image. This coiTection routine was 
applied to each reduced sky frame before they were combined to 
the master sky, and to each science frame after the subtraction of 
the corrected sky. The application of this routine after the basic 
reduction ensures that features due to dark current, detector bias 
and response are not altered by the 50 Hz noise correction, and 
are correctly removed using the dark and the twilight flat. 

Very bright sources cause a regular pattern of electronic 
ghosts, the sizes and shapes of which depend on the detector 
array and the brightness of the inducing object. The minimum 
stellar peak flux in the raw images sufficient to generate a visible 
electronic ghost was about 2700 counts for all three 2.0 s DIT 
datasets. For the dataset with the longer DIT of 20.0 s the min- 
imum flux was about 5500 counts, as ghosts of fainter stars are 
probably hidden within the brighter background. The location of 
the object on the detector aiTay determines the position of these 
ghosts. A bright star at the position (x,y) causes e lectronic ghosts 
at(x,1024-y), (1024-x, 1024-y) and (1024-x, y) (lAgeorgeset al.1 
12007 , chap. 4.7.1). For each science frame a mask covering the 
visible electronic ghosts was generated and added to the mask 
containing the bad pixels in the dark, the flat field and bad pix- 
els arising from cosmic ray hits. A few optical ghosts, appar- 
ent as sets of concentric rings with a radius o f about 40 pixel 
(=L1"), are caused by the bright stars as well ("Ageorg es et alj 
12007 , Fig. 16). Unlike the electronic ghosts, the position of these 
ghosts seems to be fixed with respect to the observed star field, 
which prohibits a coiTection for these ghosts during the image 
combination, as they are placed on the same location in the com- 
bined image. As the brightness of stars, residing within an op- 
tical ghost, cannot be reliably determined these stars were later 
removed from all source catalogues. 

The science frames of a dataset were combined into one 
combined image with the IRAF task drizzle (Fruchter& Hook 
I2002h . Each science frame was linearly weighted by the inverse 
of the full width at half maximum (FWHM) of an unsaturated, 
bright star present in all frames of this dataset in order to opti- 
mise the spatial resolution in the combined image. Pixels con- 
tained in the bad pixel masks are rejected by drizzle during 
the image combination. The large number of 44 exposures for 
the second epoch allowed for using only the 33 frames with a 
FWHM < 0.082" in order to enhance the spatial resolution with- 
out losing stars at the faint end. 

Stellar fluxes an d positions were d etermined with the 
starfinder algorithm (iDiolaiti et al.l l2000l) . which is designed 
for high precision astrometry and photometry on AO data of 
crowded fields. The point spread function (PSF) is derived em- 
pirically from the data by median superposition of selected stars 
after subtraction of the local background and normalization to 
unit flux. Using an empirical PSF is preferable for astrometric 
AO data, as the steep core and wide halo are not well reproduced 
by analytic functions. Stars whose peak fluxes exceed the linear- 
ity limit of the detector and are included in the list of stars for the 
PSF extraction are repaired by replacing the saturated core with 
a replica of the PSF, scaled to fit the non-saturated wings of the 
staiQ. Only if the saturated stars are repaired they are definitely 



' Stars, whose peak flux exceed the linearity limit of the detector are 
referred to as saturated for the remainder of this paper. 



detected and fitted by the algorithm, so that their contribution on 
the flux of neighbouring stars can be subtracted. This is of spe- 
cial importance for faint stars located within the halo of a satu- 
rated star in order to measure their fluxes precisely. As the option 
for a spatial variation of the PSF due to uncoiTected atmospheric 
modes for larger distances to the AO guide star is not yet imple- 
mented in the starfinder code, a spatially constant PSF was used. 
In order to obtain a valid estimate for the PSF for the most parts 
of the image, preferably isolated, bright stars uniformly spread 
across the image were selected for the PSF extraction. All satu- 
rated stars were included in the list of PSF stars in order to be 
repaired. The total number of selected PSF stars and the number 
of saturated stars among them are listed in Table |2] The compa- 
rably small number of saturated stars of the last dataset is due to 
the higher linearity limit of the Aladdin3 detector. 

The simplification of a constant PSF across the whole im- 
age led to spatially varying PSF fitting residuals and in turn to 
small-scale zeropoint variations across the field. This is typical 
for AO data and is mostly a consequence of anisoplanatism at 
increasingly larger distances from the natural guide star. As the 
extracted PSF resembles an average of the different PSFs across 
the image, the variation of the residuals after PSF subtraction 
is not centred at the position of the guide star. It rather shows 
a radial trend also depending on the respective image segment. 
In order to coiTect for these variations a spatially varying cor- 
rection factor was determined from the flux ratio of the residual 
flux in the PSF subtracted image and the stellar flux within an 
aperture around the centroids of isolated stars. Only stars with- 
out neighbours brighter than 1% of the star's flux were selected. 
The radius of the aperture was chosen to be either 0.08", 0.14" 
or 0.22" in dependence of the stellar flux to make sure that on 
the one hand the bulk of the residual and stellar flux is included, 
while on the other hand the residual flux for the faint stars is not 
dominated by the background noise. The flux ratio was fitted in 
dependence of the distance to the image centre for segments of 
45° either by a constant offset or a small linear trend. The coiTec- 
tion factor fconjir) for the i* segment, which is to be multiplied 
to the fluxes of all stars in this segment, follows from the fit of 
the flux ratio 'Rfujir): 



The error of fcon/f) is derived from the fitting error of 'Rnt^i(r), 
which is A'Rfitj(r) - Ac,- if the flux ratio was fitted by a constant 
off'set c/ and A^fit,^'") - ^(r^bj)^ + (Ac,)^ if the flux ratio was 
fitted by a linear trend with 'Rntj(r) = bir + a: 

^fconj(r) = f , j ■ (2) 

\l-'Rf,ij(r)l 

This procedure resulted in the most consistent photometric cali- 
bration across the observed field. 

3. Photometric calibration and error estimation 

Reference sources for the photometric cal ibration were take n 
from the Galactic Plane Survey (GPS) dLucas et alJ l2008h . 
a part of the UKIRT Infrared Deep Sky Survey (UKIDSS) 
([Lawrence et al.ll2007h . Magnitudes of stars within the UKIDSS 
catalogue are determined from aperture photometry using an 
aperture radius of I" and are calibrated usi ng the Two Micron 
AU Sky Survey (2MASS) (Sk rutskie et al.ll2006l) . Data from 
the Sixth Data Release (DR6) for the Quintuplet cluster was 
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retrieved from the UKIDSS archive (iHamblv etalJl2008h . For 
the NACO H-band data (2003) and the second epoch /Tj-band 
dataset (2008), sources from the UKIDSS catalogue in the 
Quintuplet Cluster served as zeropoint reference. For a set of 
calibration stars (29 in H-, 13 in /Tj-band), which could unam- 
biguously be assigned to calibrated sources in the UKIDSS cat- 
alogue, the individual zeropoints were determined. Due to the 
high spatial resolution of the NACO data usually several fainter 
stars can be resolved within the UKIDSS 1" aperture around the 
calibration star. As these stars do contribute to the measured flux 
in the UKIDSS survey, the PSF-flux of all stars falling within 
a radius of r = 1" - 0.5 x FWHMpsF, with FWHMpsF be- 
ing the FWHM of the extracted PSF of the respective NACO 
dataset, was added and compared to the magnitude of the re- 
spective star in the UKIDSS catalogue. The final zeropoint was 
then determined from the average of the individual zeropoints of 
the calibration stars. The zeropoints of the two /Tj-band datasets 
from the first epoch were determined subsequently using the 
calibrated second epoch data. No significant colour terms were 
found between the NACO H, K, and the UKIDSS H,K filter 
systems. 

The estimation of the photometric an d astrometric uncer tain- 
ties follows the approach described in iGhez et a D (l2008l ) and 
IStolte etalJ (l2008). The reduced science frames for each dataset 
were divided into three subsets of comparable quality and cover- 
age. Each subset of 5 (first epoch) or 1 1 frames (second epoch) 
was then combined with drizzle and the photometry and astrom- 
etry of the resulting auxiliary image was derived with starfinder 
in the same way as for the deep images. Each auxiliary frame 
was calibrated with respect to the deep image. As we are only 
concerned about the three independent measurements no correc- 
tion factor fcomjir) for the flux was applied, such that the pho- 
tometric uncertainty derived from the auxiliary frames includes 
only the PSF fitting uncertainty. The photometric and astromet- 
ric uncertainty was derived as the standard deviation of the three 
independent measurements for each star detected in all three 
auxiliary frames. As no preferential direction is expected for the 
positional uncertainty, the astrometric uncertainty of each star 
is computed as the mean of the positional uncertainty in the x- 
and y-direction. The astrometric and photometric uncertainties 
as derived from the auxiliary frames are shown in dependence of 
the magnitude in Fig. |2] for all datasets. 

In order to remove false detections from the three /Tj-band 
catalogues, only stars which were detected in all three auxiliary 
images of a dataset, and hence with measured astrometric and 
photometric uncertainties assigned to them, were kept in the re- 
spective source catalogue. For the //-band data this criterion was 
not applied. The //-band was matched in the further analysis (see 
Sect.|6]l with a /Tj-band catalogue containing only stars detected 
in both epochs. It is assumed that a star found in the /Tj-band 
images of both epochs is a real source and if it is missing in 
one of the //-band auxiliary images this is a consequence of the 
substantially lower photometric depth of the auxiliary image. 

The photometric errors as stated in the final source catalogue 
(Table HI do include the respective zeropoint uncertainties, the 
photometric uncertainties due to the flux measurement from PSF 
fitting, and the error of the correction factors as given in equa- 
tion ©. 

4. Completeness 

In order to quantify the detection losses due to crowding effects, 
the local completeness for each dataset was determined from the 
recovery fraction of artificial stars inserted into each combined 



Table 2. Number of stars for PSF extraction. 



Dataset No. No. of PSF stars No. of saturated PSF stars 

1 37 17 

2 239 136 

3 48 29 

4 69 15 



image. The artificial star experiment for the //-band data cov- 
ers a magnitude range from 9.5 to 2 1 .5 mag. For each magnitude 
bin with a width of 0.5 mag, 42 artificial star fields were gener- 
ated. Each artificial star field was created by adding 100 artificial 
stars, which are scaled replica of the empirical PSF, inserted at 
random positions and with random fluxes within the respective 
flux interval, into the combined image. 

For the three /T^-band datasets, the artificial stars were in- 
serted at the same physical positions as in //-band and with a 
magnitude in Kg yielding a colour for the respective artificial 
star of H - Ks = 1 .6 mag, which resembles the colour of main 
sequence (MS) stars in the Quintuplet cluster (see Sect. |6]l. The 
photometry on the images with added artificial stars was per- 
formed in the same way and using the same PSF as for the 
original images. In addition to artificial stars which were not 
re-detected by starfinder, also stars where the recovered magni- 
tudes deviated strongly from the inserted magnitudes, were con- 
sidered as not recovered. The criterion to reject recovered stars 
due to their magnitude difference between input and output mag- 
nitude was derived from polynomial fits to the median and the 
standard deviation of the magnitude difference within magnitude 
bins of 1 mag (Fig. 13). Extreme outliers with absolute magnitude 
differences larger than 20 times the standard deviation were ex- 
cluded from the determination of the median and standard de- 
viation used for the fits. Stars with absolute magnitude differ- 
ences exceeding 0.20 mag and being larger than 1 .5 times the fit 
to the standard deviation are treated as not recovered. The me- 
dian of the magnitude difference exposes a systematic increase 
towards the faint end, exceeding 0.05 mag for Kg > 19.4 mag or 
H > 20.25 mag. This trend indicates that for the faintest stars 
the measured fluxes are not reliable anymore. As we restrict the 
analysis to stars brighter than Ks < 19 mag, sources at these faint 
magnitudes are excluded from the proper motion and mass func- 
tion derivation. 

The left panel in Fig.|4]shows the overall recovery fraction in 
the /Tj-band for the two epochs and the //-band within a radius of 
500 pixels (=13.6") from the image centre, the part of the image 
actually used for the determination of the present-day mass func- 
tion (see Sect.O. The shown recovery fraction for the /T^-band 
data from 2003 is a combination of the recovery fractions for the 
two /Tj-band datasets of that epoch. The dataset with the longer 
DIT of 20.0 s is used only for magnitudes fainter than the linear- 
ity limit of this dataset at 14.3 mag. For brighter magnitudes, the 
recovery fraction of the 2003 /T^-band data with the short DIT 
of 2.0 s is drawn. The star catalogues of these two datasets are 
combined accordingly in the subsequent analysis (see Sect. 15.2b . 
The total recovery fraction also shown in the figure is the product 
of all three recovery fractions and is most relevant for the com- 
pleteness correction of the mass function, as i) only stars which 
are detected in both epochs can be proper motion members and 
ii) only for stars with measured //-band magnitudes can masses 
be derived. 

Completeness varies as a function of position due to the non- 
uniform distribution of brighter stars and hence is a function 
of the stellar density and magnitude contrast between neigh- 
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Table 1. Overview of the VLT/NAOS-CONICA datasets. 



Dataset No. 


Date 


Filter 


No. of Frames 


DIT 


NDIT 


Lint 


Airmass 


Seeing 


FWHM of PSF 


Strehl ratio 










(s) 




(s) 




(") 


(") 




1 


2003-07-22 


H 


16 


2.0 


30 


960 


1.00-1.02 


0.47-0.60 


0.078 


0.15 


2 


2003-07-22 




16 


20.0 


2 


640 


1.03-1.06 


0.36-0.49 


0.080 


0.22 


3 


2003-07-23 




16 


2.0 


30 


960 


1.03-1.07 


0.31-0.47 


0.082 


0.26 


4 


2008-07-24 




33 


2.0 


15 


990 


1.00-1.01 


0.49-0.60 


0.080 


0.26 



Notes. Total integration time of the central part of the image with maximum overlap. 



hours (see e.g., lEisenhauer et al.l 4998': 'G ennaro et al.ll20 l7). A 
spatially-dependent approach to determine the local complete- 
ness value becomes especially important if the cluster exhibits a 
non-symmetric geometry or in the presence of sparse very bright 
objects, which heavily affect the completeness values in their 
surrounding as in the case of the Quintuplet cluster In order 
to assign a local completeness valu e to each detec t ed sta r, the 
method described in Appendix A of iGennaro et al. was 
applied to derive completeness maps for each combined image 
containing the recovery fraction for every pixel as a function of 
magnitude. The procedure encompasses threes steps performed 
for each m agnitude bin (for a d etailed description the reader is 
referred to iGennaro et al.ll201 ih : 1.) Derive the local, averaged 
completeness value at the position of each artificial star from the 
number of recovered nearest neighbours including the star itself. 
2.) Interpolate these local completeness values into the regular 
grid of image pixels. 3.) Smooth the obtained map with a box- 
car kernel with a width of the sampling size in order to remove 
potential artificial features introduced by the previous step. For 
the last step the completeness maps of all magnitude bins are 
used. To ensure that the completeness decreases monotonically 
with increasing magnitude, a Fermi-like function is fitted to the 
completeness values at every pixel in the image as a function 
of magnitude. The completeness (or recovery fraction) for every 
real star in the respective band and dataset can then be com- 
puted from the fit parameters at the position of the star. The right 
panel in Fig. |4] shows the combined Ks-hand image for the sec- 
ond epoch with superimposed 50%-completeness contours stat- 
ing the limiting magnitudes. The very bright stars with their ex- 
tended halos hamper the detection of nearby faint stars causing 
the recovery fraction to be non-uniform across the field, as ex- 
pected. The completeness of a star entering the mass function is 
the product of its completeness in the //-band, the second epoch 
/Tj-band data and either in the 2.0 s DIT (/Tj ^oos < 14.3 mag) or 
the 20.0 s DIT (/T, 2003 > 14.3 mag) first epoch /Tj-band dataset 
as determined from the respective completeness maps: 

/comp ~ /comp,Ks2{)08 ^ ,/comp,K,s2003 ^ ./comp, H2003 ■ (3) 

For stars brighter than // = 1 3 .5 mag 01 Kg - 1 0.4 mag the com- 
pleteness was assumed to be 100%. 

5. Proper motion membership 

Due to the high field star density for lines of sight towards the 
Galactic centre the distinction between cluster and field stars 
becomes particularly important. As most of the field stars are 
located within the Galactic bulge they have similar extinction 
values as the cluster and cannot be distinguished from cluster 
members on the basis of their colours alone. The high astromet- 
ric accuracy of the AO assisted VLT observations in combina- 
tion with the time baseline of 5.0 yr allows for the measurement 
of the individual proper motions of stars at the distance of the 



Quintuplet cluster. The primary applied method to discern the 
cluster members from the field stars is based on the measured 
proper motions. 

5.1. Geometric transformation 

In order to determine the spatial displacements, two geomet- 
ric transformations were derived to map each position in the 
two first epoch Kg images (2003) with short (2.0 s) and long 
(20.0 s) DIT onto the correct position in the second epoch K, 
image (2008). The second epoch is used as reference epoch be- 
cause of the higher astrometric accuracy, deeper photometry and 
brighter linearity limit of this dataset. Only the /T^-band datasets 
were used to determine the spatial displacements, as due to their 
higher Strehl ratios the stellar cores are better resolved than in H- 
band, providing the better centroiding accuracy and hence most 
accurate astrometry. 

Under the assumption that internal motions are not resolved 
so far, the cluster itself served as the reference frame. The 
geometric transformation was derived in an iterative process. 
First, a rough transformation was determined with the IRAF 
task geomap using the positions of manually selected bright, 
non-saturated stars uniformly distributed across the images of 
both epochs. The respective catalogue of the first epoch dataset 
(2003) was then mapped onto the catalogue of the second epoch 
(2008) to get a mutual assignment of stars found in both cata- 
logues. From these stars the most likely cluster candidates were 
selected to provide the reference positions for the refined, fi- 
nal geometric transformation. As the bright stars used for the 
preliminary transformation are likely cluster members the dis- 
tribution of spatial displacements in the x-,y-direction of clus- 
ter star candidates are expected to scatter around the origin. 
Therefore only stars with spatial displacements within a radius 
of 0.8 pixel=4.3 mas/yr from the origin were selected for the 
derivation of the final transformation, which excludes most of 
the presumed field stars. Further, as the bulk of cluster stars are 
probably brighter than most stars in the field, only non-saturated 
bright and intermediate bright stars (11 .5 < Ks < 15.5 mag for a 
DIT of 2.0 s and 14.0 < K, < 17.0 mag for a DIT of 20.0 s) pro- 
vide the reference positions. The final geometric transformations 
were derived with geomap in an interactive way. The residual 
displacements in the x-, y-directions between the transformed 
first epoch and the second epoch coordinates were minimized 
by iteratively removing outliers and carefully adapting the order 
of the polynomial fit (= 3 for the final transformations). The final 
rms deviation of the geometric transformation was 0.2 mas/yr in 
the X- and 0.3 mas / yr in the y-direction for the dataset with a DIT 
of 2.0 s, and 0.3 mas/yr in the x- and y-direction for the dataset 
with a DIT of 20.0 s. For a total cluster mass of Md ~ 5900 M© 
within a radius of r < 0.5 pc (see Sect. |8]l the internal velocity 
dispersion is expected to be on the order of 0.15 - 0.2 mas/yr or 
6 - 8km/s. As this is smaller than the uncertainty of the geo- 
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metric transformation alone, intrinsic motions are not resolved. 
Therefore the above assumption is justified. 



5.2. Data selection and combination 



Each of the two transformed star catalogues of the /Tj-band 
data from the first epoch was matched with the star cata- 
logue of the second epoch using a matching radius of 4 pixels 
(=108 mas - 1.4 x FWHMpsp)- The matching radius was cho- 
sen small enough to avoid mismatches between close neigh- 
bouring stars, but large enough to include all moving sources 
at GC distances below the escape velocity of the GC. A dis- 
placement of 108 mas within the time baseline of 5.0 yr for a 
distance of 8.0 kpc to the GC dGhez et al.ll2008l) corresponds to 
a proper motion of 820km/s. The combined astrometric uncer- 
tainty tTpo, = ^'^pos,Ks2003 + '^pos.Ksaoos ^as derived for both of 
these catalogues. In the left panel of Fig. |5] the combined as- 
trometric uncertainty is plotted against the magnitude. The dat- 
apoints below the linearity limit of the long exposure in 2003 at 
Ks - 14.3 mag originate from the match of the second epoch 
with the data with a DIT of 20.0 s, the datapoints at brighter K, 
magnitudes are from the match with the data obtained with a 
shorter DIT of 2.0 s. For the matched catalogue using the first 
epoch dataset with a DIT of 20.0 s, the median and the stan- 
dard deviation of the astrometric uncertainties within bins of 
0.5 mag width were fitted by a third and second order polyno- 
mial, respectively. The fit to the median and the sum of both 
fits are shown in all three panels of Fig. |5] for comparison. The 
usage of one averaged PSF for the whole image results in the ob- 
served radial increase of the PSF fitting residuals (see Sect. |2]i. 
The centroiding accuracy is therefore expected to decrease to- 
wards larger radii resulting in a larger astrometric uncertainty. 
The centre and right panel of Fig. |5] exemplify this behaviour 
by using only stars with a distance of less than or greater than 
500 pixel = 13.6" from the centre of the combined images in 
both epochs, respectively. The decrease in the scatter and mag- 
nitude of the astrometric uncertainties is striking. The median of 
the astrometric uncertainty for 12 < /Tj < 18 mag is 2.52 mas for 
stars within a radius of 500 pixel, but 4.58 mas for stars outside 
that radius. Therefore, the further analysis is restricted to stars 
within a radius of 13.6" from the centre of the observed field 
of view for the remainder of this paper. The astrometric uncer- 
tainties rise steeply near the detection limit at about 20 mag (see 
centre panel in Fig. |5]l. For stars fainter than 19 mag, almost no 
stars exhibit an uncertainty below the median value of stars with 
intermediate brightness (14 < Kg < 17 mag). Stars with a Kg- 
band magnitude fainter than 19 mag are therefore excluded from 
the sample. As last selection based on the combined uncertainty, 
stars fainter than Ks - 14.3 mag are removed if their uncertainty 
is above the sum of the fits of the median and standard devia- 
tion derived from the combined uncertainty of all observed stars 
(see Fig. |5]l. The percentage of rejected stars varies between 
and 9.5% for the affected magnitude bins and does not show a 
systematic trend with magnitude, therefore no systematic bias is 
introduced by this selection. After the above mentioned selec- 
tions, the two matched catalogues were combined. Stars fainter 
than Ks - 14.3 mag were taken from the match with the DIT 
20.0 s first epoch data, brighter stars originate from the matched 
catalogue using the dataset from the first epoch with a DIT of 
2.0 s. The final /Tj-band catalogue contains a total of 1297 stars. 



5.3. Ttie proper motion diagram 

Individual stars are plotted in the proper motion diagram (Fig.|6l) 
with proper motions in the east-west-direction on the x- and 
proper motions in the north-south direction on the y-axis. As the 
cluster is used as the reference frame, the distribution of clus- 
ter members is centred around the origin and overlaps with the 
elongated distribution of the field stars. The orientation of the 
field stars is approximately parallel to the plane of the Galaxy 
(dashed line in Fig.|6]l. The dotted line running through the ori- 
gin and vertically to the Galactic plane splits the proper motion 
diagram into two halfs being referred to as the North-East seg- 
ment (upper half) and the South- West segment (lower half). 

Figure Q shows histogram plots of the distribution of the 
proper motions parallel (left panel) and vertical to the Galactic 
plane (centre panel). The distribution of proper motions in the 
direction parallel to the Galactic plane is strongly peaked at the 
origin, with a very steep decline in the North- West segment, and 
a slightly broadened decline and overlap with the broad field star 
distribution in the South- West segment, as expected from Fig.|6] 
The proper motions vertical to the Galactic plane are almost dis- 
tributed symmetrically with respect to the Galactic plane (centre 
panel in Fig.|7]), confirming the assumed orientation of the field 
star distribution in the proper motion diagram. This and the ex- 
posed offset of the field star distribution in the proper motion 
diagram indicate a movement of the Quintuplet Cluster paral- 
lel to the Galactic plane towards North-East with respect to the 
field as was found previously for the Arches Cluster ( Stolte et aH 
2008). The sample of stars with proper motions in the North- 
East segment is least contaminated by field stars and was there- 
fore used to derive the membership criterion. The distribution 
of proper motions in the North-East segment was fitted with a 
Gaussian function (right panel in Fig. |7]i. Stars whose proper 
motions are within a circle of radius Icr - 2.24mas/yr, where 
cr is the width of the Gaussian fit, are selected as cluster mem- 
bers (see Fig.lHl. Two of the initial five Qui ntuplet members (Ql, 
Q9) (Na gata et al.ll990l:IOkuda et aljl 990) do not fall inside this 
circle. Their fluxes are exceeding the linearity limits by a factor 
of 8-30, such that their positions are not well determined. Note 
that this only affects the very brightest sourc es, for which spec- 
troscopic mem ber identification is available jFiger et al]|1999bt 
iLiermann et a l. 2009). These two stars were added manually to 
the sample of proper motion members. 

6. Colour-magnitude diagrams 

The Ks source catalogues of proper motion members and non- 
members were matched with the source catalogue of the first 
epoch H-band data. All 1218 matched stars (member and 
non-members) are included in the final source catalogue (see 
TablelU. The corresponding colour-magnitude diagrams (CMD) 
are shown in Fig. |8] and use only magnitudes in H and Ks from 
the first epoch to avoid additional scatter being introduced by 
variable stars. Stars whose fluxes exceed the respective linear- 
ity limit in either H and/or Ks at H = 12.05 mag and Kg = 
1 1.25 mag are marked with crosses. All 1218 stars (member and 
non-members) with measured proper motion and (H-Ks) colour 
are included in the final source catalogue (see Table|4|. 

Cluster and field stars separate well as can be seen by char- 
acteristic features of the field population (right panel), that are 
absent in the cluster selection (left panel). For example an elon- 
gated overdensity is observed, which starts at about H - 17 mag 
, H - Ks - 1.8 mag and extends to redder c olours along the 
reddening path adopting the extinction law by iNishivama et al.l 



6 



B. HuBmann et al.: The present-day mass function of the Quintuplet cluster 



(I2009h . It is consistent with arising from red clump stars located 
in the Galactic bulge. Assuming the intrinsic ^T-band magnitude 
for red clump stars of A' - -1.61 mag by Alves (2000), the as- 
sumed intrinsic colour of H — Ks - 0.07 of Nishiva ma et aP 
(|2006|) . a distance to the Galactic Centre of 8 kpc (iGhez et al.l 
|2008|) and an approximate extinction for the cluster of Ak^ - 
2.35 mag yields Hrc - 17.05 mag and (H - Ks)rc = 1-79 mag. 

Several blue foreground stars with colours H-K, < 1.3 mag 
are seen to the left of the cluster member sequence (Fig. [8] left 
panel). These sources are likely disc main sequence stars fol- 
lowing the differential rotation of the outer Milky Way rota- 
tion curve. With expected velocities of ~ 200km/s, they can- 
not be distinguished from the cluster population on the basis 
of their proper motion alone. Furthermore a few very red ob- 
jects, which could be non-members by comparison with the 
field CMD, remain in the proper motion sample. In order to re- 
move these contaminants a two-step colour-cut was applied to 
stars fainter than H - 14 mag. First the blue foreground and 
red background stars were removed by keeping only stars with 
1.3 < H - Ks < 2.3 mag. In a second step the individual ex- 
tinction of the remaining stars fainter than H = 14 mag was de- 
termined from the intersections of the lines of reddening with a 
4 Myr isochrone assuming a distance to the cluster of 8 kpc. The 
method to derive the individual extinction and the used isochrone 
are explained in detail in Sect. |7] The isochrone was shifted to 
an extinction of A^^ = 2.88 mag, corresponding to the sum of 
the mean (Ak^ = 2.41 mag) and twice the standard deviation 
(cT/i^ = 0.24 mag) of the individual extinctions of the cluster 
members remaining after the first colour-cut, and stars redder 
than the shifted isochrone were also removed from the sample 
of cluster stars. 

The designated cluster members and non -members were 
compa red with the K-band spectral catalogue of lLiermann et alj 
(l2009h (further abbreviated as LHO catalogue) for a spectral 
classification of the brighter stars and in order to assess the selec- 
tion of cluster stars based on their proper motions and colours. 
Only observed stars with a Kg band magnitude brighter than 
15.5 mag, which is about 1 mag fainter than the faintest star in 
the LHO catalogue, were included in the comparison. Eighty- 
five stars from the spectral catalogue could be assigned to 92 
observed stars (69 members, 23 field stars). The ambiguous as- 
signments of 6 stars from the LHO catalogue to 13 observed 
stars (all members) are caused by the lower spatial resolution of 
the SINFONI-SPIFFI instrument of 0.250" for the used 8" x 8" 
field of view. The spectral classification for the matched stars 
is indicated in a simplified form by the overplotted symbols in 
Fig. [8] Fig. |9] and Fig. [TO] Stars with ambiguous assignments 
are additionally mar ked with an X-c r oss. O ne star (LHO 110) 
was re-classified in iLiermann et a D (1201 Oh from 06-8 I f to 
WN9h and is treated accordingly in the figures. The numbers 
and spectral classifications from the LHO catalogue are noted in 
the source catalogue (Table|4]i. Six late-type M,K-giants are still 
contained within the cluster sample after the colour-cut and are 
very likely remaining contaminants with motions similar to the 
cluster members from the Galactic bulge considering the young 
age of the cluster These stars and stars rejected by the colour- 
cut were removed from the final cluster sample and added to the 
proper motion non-members in the field star CMD (plotted as 
triangles in the right panel of Fig. [8]l. The one early-type star 
(04-7 I f) among the designated field stars is located at the edge 
of the analysed area of the data and very close to a second star 
just outside this region. It is therefore unclear if the spectral clas- 
sification really belongs to this star or its neighbour, hence the 
star was not added to the final cluster sample. For 12 of the 62 



stars in the LHO catalogue, which could be assigned to desig- 
nated cluster members, the sources in our catalogue exceed the 
/Tj-band linearity limit by more than 1 mag. This impedes the 
repair of the core by starfinder and the correct measurement of 
the position and proper motion. Disregarding these 12 stars, the 
percentage of contaminating M,K giants, which cannot be dis- 
cerned from the cluster members based on their proper motion 
or colour, amounts to 6/(62 - 12) = 12%. 

The CMD of the final sample of cluster stars is shown in 
Fig. |9] and, separated into the North-East and South- West seg- 
ment of the proper motion diagram, in Fig.[TOl The slight over- 
density located at H - - 1.8 mag, H - 17 mag indicates a 
remaining contamination with red clump stars, which is more 
pronounced for stars with proper motion in the South- West seg- 
ment. The CMD for the South- West segment contains 92 stars 
more than for North-East segment mainly at the faint end of the 
observed population, which appears slightly broadened. This is 
expected from the proper motion diagram as the field star popu- 
lation overlaps with the cluster stars in the South- West segment 
causing a larger contamination for this segment. The astrometric 
uncertainty and therefore the scatter in the proper motion dia- 
gram increases for fainter magnitudes and therefore the confu- 
sion with faint field stars is more severe. The cluster members in 
the North-East segment therefore constitute the cleanest sample. 



7. Mass derivation 

Based on the presence of WC stars, O I stars and a red su- 
pergiant within the Quintuplet cluster, iFiger et aTl (Il999bl) de- 
rived an average age of 4 + 1 Myr assuming a coeval popula- 
tion. More recently the ages of 5 WN stars were determined 
by comparison of their luminosities and effective temperatures 
as derived from spectral line fitting with stellar evolution mod- 
els to be about 2. 4 - 3.6 Myr pointing to a somewhat younger 
age of the cluster (ILiermann et aLllToiOh . To study the influence 
of the assumed cluster age on the slope of the mass function, 
three isochrones with ages of 3, 4 and 5 Myr were used to de- 
rived the initial stellar masses. The isochrones are a combina- 
tion of Padova main sequence (MS) isochrones and pre-main 
sequence (PMS) iso chrones derived froni Pisa-FRANEC PMS 
stel lar models (see 'Genna ro et"ai] 1201 It iMarigo et all 120081: 
DeglTnnocenti et al. 2008). As the NACO photometry is cah- 
brated by means of UKIDSS sources (see Sect.O, the combined 
isochrones, for simplicity referred to as 3, 4 and 5 Myr Padova 
isochrones in the following, were transformed from the 2MASS 
into the UKID SS photometric system using the colour equations 
(6) - (8) from iHodgkin et all (12009) . To cover the effect of a 
different set of stellar models on the derived masses, a 4 Myr 
Geneva MS isochrone with enhanced mass loss for high mass 
stars, M > 15 Mq, (Leieune & Schaerer 2001) was included 
in the comparison. The conversion of this isochrone into the 
UKIDSS filter system encompa ssed two steps. The iso chrone 
was first transformed from the iBessell & Brett! ( Il988h to the 
2MASS photometric system using the updatecfl transformation 
by Carpenter (2001) and subsequently from the 2MASS to the 
UKIDSS filter system using the above mentioned conversion. 

For all isochrones, solar metallicity according to the de- 
scription of the underlying stellar modelfl was assumed, and 



- Carpenter, J.M., 2003 see^ttp: //www. astro. caltech. edu/~ jmc/2mas; 
^ solar metallicity for the G eneva isochrone: X = 0.68, Y = 0.3, 
Z = 0.020 jLeieune & Schaeredl2001.) : 

solar metallicity for the Pad ova isochrones: X = 0.708, Y = 0.273 , 
Z = 0.019 for M < IMq dMarigo et al.ll2008l : iGirardi et alJl2000h . 
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a distance to the Galactic centre of 8.0 kpc (iGhez et alJl2008l) 5 
was applied as the distance to the Quintuplet cluster. The 
four isochrones shown in Fig. [TT] were reddened by a fore- <' 
ground extinction of A / c = 2.35 mag using the extinction 
law of iNishivama et all (f2009l) (Ah : Ak^ = 1.73 : 1) to | 
match the observed MS of the cluster members. This extinc- o 
tion law is one of the most recent determinations of the ex- " 
tinction in the near-infrared along the line of sight towards the 
G alactic centre and consistent wi th other current fin dings, e.g., 
bv lStraizvs & Laugalvsl (l2008h or lSchodel et"an (1201 Q) . 

The individual mass and extinction of each star in the final 
cluster sample was determined from the intersection of the line 
of reddening through the star with the respective isochrone in the 
CMD. Due to the local maximum of the PMS at the low-mass 
end as well as the extended loop at the transition from the end of 
the hydrogen core burning to the contraction phase at the high- 
mass end, the de-reddening path of a star may have several inter- 
sections with the isochrone, thus leading to an ambiguous mass 
assignment (the affected areas in the colour-magnitude plane are 
shaded in grey in Figs. l9l[T0l and fTTTi. For these stars, the masses 
at each intersection were averaged. The post-MS phase after 
the exhaustion of hydrogen in the stellar core is very rapid (a 
few 10^ yr according to the stellar models) and apparent in the 
isochrones as the branch with increasing //-band brightnesses, 
re-rising after the decline connected to the contraction phase. 
Due to its short duration, which causes the Hertzsprung gap in 
the Hertzsprung-Russel diagrams of stellar clusters, only the two 
intersection points with the upper part of the MS and with the 
subsequent falling branch of the isochrone were averaged. Two 
O stars from the LHO catalogue have no intersection with the 
Geneva isochrone on the MS or the falling branch, therefore an 
initial mass of 47.3 M0, which is the maximum mass along this 
isochrone used for the mass determination (see Table |3]l, was 
assigned to them. 

1 1 Wolf-Rayet stars out of the 21 observed in the Quintuplet 
Cluster are contained within our sample of cluster members. 
The masses for these stars could not be determined from the 
isochrones but the mass range of Wolf-R ayet stars was i nferred 
from the underlying stellar models by |B ressan et al.l dl993h 
for the Padova iso chrones and bv iMev net et al. (1995) and 
ISchaller et al.l (1 19921) for the Geneva isochrone (see TableO. 

Figure [12] shows an enlarged detail of the transition re- 
gion from the PMS to the MS in the member CMD. The three 
Padova isochrones shown are all shifted to an extinction of 
Ak^ - 2.35 mag in order to fit the cluster MS. Although stars 
scatter across the whole transition region, there is an indication 
of a slight accumulation of stars at the turn-over of the PMS for 
a cluster age of 4 Myr. At this stage of evolution the CN cycle 
supports most of the stellar radiative losses, but the NO cycle is 
still not in complete equilibrium. After suflicient N is produced, 
the star undergoes a rapid contraction and reaches the MS after 
the CN-NO cycle is in complete equilibrium. The MS-turn-on 
point for this cluster age is marked by a second overdensity of 
stars (see Fig.fT2ll. In contrast, the turn-over of the PMS and MS- 
turn-on point of the 5 Myr isochrone as well as the MS-turn-on 
for the 3 Myr isochrone are located in a region almost devoid 
of stars. Irrespective of the definition of the "true" age of the 
cluster, we therefore consider the 4 Myr Pisa-FRANEC/Padova 
isochrone as the best fitting isochrone within the investigated set 
of models. 
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Table 3. Summary of isochrone properties relevant for the mass derivation. 



Isochrone name Description Wolf-Rayet mass 



3 Myr Padova 

4 Myr Padova 

5 Myr Padova 
4 Myr Geneva 



MS" + PMS* 
MS" + PMS* 
MS" + PMS* 

MS'- 



(Mq) 
85 - 100 
51-65 
37-40 
48-60 



stars with ambiguous mass assignments 



PMS 
mass range 
(Mq) 

2.1-5.5 
1.9-4.6 
1.8-4.0 



-> MS 

No. of stars 

145 
145 
149 



MS -» 
mass range 

(Mq) 
42.3 - 84.4 
39.9-50.8 
33.5 - 35.4 
45.8 -47.3 



post-MS 
No. of stars 

11 

5 
4 




Notes. Padova isochrone with solar metallicity for m > 4M0 jMarigo eL alJl200§) . Pre-main sequenc e parts of the isochrones (m < 4Mq) 
are derived from Pisa-FRANEC PMS stellar models (Degl'Innocenti et al. ( 20OT). see lGennaro et al.l ( 120111) for t he combination with the Pa dova 
isochrones). Geneva isochrone with solar metallicity and enhanced mass loss for high mass stars, M > 15 Mq, jLeieune & Schaere3l200lh . 



8. Mass functions 

In order to avoid potential biases introduced by bins with a 
very small number of objects or large differences in the num- 
ber of stars between the low- and high-mass bins, we adopted 
the method proposed bv lMafz Apellaniz & Ubedal (l2005l) . Here, 
the widths of the different bins are adjusted such that each bin 
houses approximately the same number of stars (Method A). If 
the number of stars in the sample did not split up evenly for the 
chosen number of bins, the bins to contain one additional star 
from the remaining stars were chosen randomly. The stars were 
then sorted according to their masses and distributed among the 
bins. For each isochrone the mass function (MF) and slope were 
determined for dividing the cluster sample into 4, 8, 12, 16 and 
20 bins. The boundary between two adjacent bins was set to the 
mean of the most/least massive star in the respective bins. The 
minimum mass used for each mass function was set to the low- 
est mass of a star with a unique mass assignment for the respec- 
tive isochrone, i.e. lying above the ambiguity region caused by 
the PMS/MS transition. Stars with ambiguous mass assignments 
at the upper end of the MS were kept for the mass function, 
as due to their small number they all contribute to the upper- 
most bin in the mass function. The upper mass limit or upper- 
most bin boundary m^p was calculated from the data to be (see 
iMafz Apellaniz & Ubedal (l2005h ) 

mup = m„ -H 0.5 (m„ - m„_i) , (4) 

with n being the total number of stars. The number of stars in 
each bin «, was normalized by the respective bin width Am, . The 
logarithm of the normalized number of stars per bin as a function 
of the logarithm of the medium mass of each bin was fitted with 
a straight line using the IDL routine UNFIT, which performs a 
minimisation. 

The uncertainty of the number of stars per bin An, is derived 
bv Mafz Apellani z & Ubedal (l2005h from the standard error of a 
binomial distribution («p, (l - p,))''^, where the unknown true 
probability for a star to reside in the /* bin p, is approximated 
by the measured value 

An,^J'l±Z^. (5) 
\ n 

Note that this uncertainty differs from the Poisson error y/rTj, 
which is usually applied to binned data. For the linear fit each 
bin was weighted by its statistical weight w,- = 1/A«,-^. The 
statistical weight w, assigned to the logarithm of the normal- 
ized number of stars per bin (log[o(n,7Am,)) follows from error 



propagation of An,- in the logarithmic plane (see equation (7) in 
iMafz Apellaniz & Ubedall2005l) : 



n - Hi 

It is basically the same for every bin as the number of stars per 
bin varies by a maximum of one. 

Besides the binning method just described the mass function 
was also determined using an equal logarithmic width for each 
bin (Method B), which is still the most common binning method 
for der iving mass function slopes (see lMafz Apellaniz & Ubedal 
( |2005|) for a discussion of the biases of this method). The lower 
and upper mass limits were determined in exactly the same way 
as above and the logarithmic bin widths were set by dividing the 
so defined mass range into 4, 8, 12, 16 and 20 bins. The weights 
applied to each bin were again calculated with equation (|6]l and 
are decreasing going to higher masses due to the lower number 
of stars contained in the high mass bins. In order to study the 
influence of the weights on the slope for this binning method, 
the slope was derived from a linear fit to the MF with and 
without weighting. 

The reported slopes of the mass function a refer to a power- 
law distribution in linear units (dn/dm oc m") with the st andar d 
Salpeter slope being a - -2.35 in this notation (Salpeter 19551) . 
If not mentioned otherwise, the mass function and its slope were 
determined using all cluster members from both the North-East 
and South- West segment (see Fig.fTOt and distributing the stars 
into bins with (almost) constant number of stars (Method A). All 
shown linear fits to the respective mass functions were derived 
from the completeness corrected mass function using for each 
star its individual completeness correction (see Sect.|4|i. 

The minimum mass of a star with unique mass assign- 
ment was 5.5, 4.6 and 4.0 M© for the 3, 4 and 5 Myr Padova 
isochrone, respectively. The minimum mass for the 4 Myr 
Geneva isochrone was set to 4.5 Mq in order to use exactly 
the same stars as for the Padova isochrone of the same age. As 
mentioned in Sect. |7] it was not possible to infer the individ- 
ual masses of the WR stars from the isochrones. Therefore, a 
constant mass within the mass ranges of the Wolf-Rayet stars 
deduced from the stellar models (see Table [3]) was assigned to 
each identified Wolf-Rayet star in the cluster sample in depen- 
dence of the assumed cluster age. The uppermost bin boundary, 
calculated with equation (|4|i, is then identical to the assigned WR 
mass. The chosen WR mass has a significant impact on the de- 
rived slopes due to the fairly large mass range of the Wolf-Rayet 
stars for cluster ages of 3 and 4 Myr. A larger assigned WR mass 
biases the mass function to a steeper slope due to the normaliza- 
tion of n, by the binwidth. The maximum difference between the 
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slopes using the minimum and maximum WR masses for each 
of the isochrones was 0.23, which is about twice the typical for- 
mal fitting error of the slope. To avoid the described bias, the 
Wolf-Rayet stars were not included in the mass function. After 
the exclusion of the Wolf-Rayet stars, the uppermost bin bound- 
ary is determined by the two most massive stars in the respective 
sample (see equation]!]). 

In order to quantify the effect of the random selection of bins 
to contain one additional star from the remainder of the division 
of the total number of stars by the number of bins (Method A), 
the distribution process and the fit to the resulting MF was re- 
peated 1000 times. The reported slopes for this binning method 
are the mean slope of all these repetitions. The maximum differ- 
ence between the slopes of the same MF due to different random 
distributions of the surplus stars was 0.02, which is very small 
compared to the formal fitting errors. 

The slope of the mass function of each isochrone was deter- 
mined using 4, 8, 12, 16 and 20 bins. Using only 4 bins results in 
slopes being systematically shallower than for the other numbers 
of bins by up to 0.12. The choice of 20 bins introduces a bias in 
the case of the 3 Myr isochrone. Due to the large number of mas- 
sive stars with averaged masses for this isochrone (see Table ]3]), 
these stars fill up the uppermost bin completely. The mass range 
of stars with averaged masses is compressed, which in turn leads 
to a decreased binwidth of the last bin. As the number of stars 
is normalized by the bin width, the normalized number of stars 
in the last bin is increased leading to a flatter slope. The most 
rehable mass function slopes are therefore obtained using 8, 12 
or 16 bins. The maximum difference of the obtained slopes for 
a given isochrone between these three bin numbers was 0.03. 
Given this negligible influence of the bin number, all results pre- 
sented in the following are determined using 12 bins (see also 
Table]5l). 

Figure 1T3] shows the comparison of the mass functions and 
the derived slopes for distributing the data in bins with (almost) 
constant number of stars (Method A, left panel) and for us- 
ing bins of equal logarithmic width of AlogjgW = 0.084 dex 
(Method B, right panel) for the 4 Myr Padova isochrone. In the 
shown example and in general the three slopes of the weighted 
fit to the mass function derived with a uniform number of stars 
per bin and the weighted and unweighted fit to the mass func- 
tion with an equal logarithmic bin widths agree well within the 
errors, if the full sample of stars is fitted. The use of a con- 
stant logarithmic bin size with applied weights following the 
prescription of iMaiz Apellaniz & Ube da (2005) generates con- 
sistent results compared to the use of bins with variable widths 
and equal numbers of stars also if only stars from the North- 
East or South- West segment of the proper motion diagram are 
included in the mass function. In contrast the unweighted fit re- 
sponds much more sensitively to fluctuations of the number of 
stars in the higher mass bins, especially if the last bin included in 
the fit is depleted, which is the case for the South- West segment. 
For the remainder of this paper only the results determined from 
mass functions with an equal number of star per bin are consid- 
ered. 

The first row in Fig. 114] shows the mass function for the 
4 Myr Padova isochrone for all cluster members (left panel), 
for stars in the North-East-segment (centre panel), and for stars 
in the South- West segment (right panel). All three slopes agree 
well, although a lack of stars in the mass function of the North- 
East segment in the mass range of 11 - 22 M0 compared to 
the South- West segment is evident. This surplus of stars for the 
South- West segment is also apparent in the CMD (Fig. [TOl i and 
could indicate a remaining contamination with red clump stars. 



Further contaminations suggested by the difference in the num- 
ber of cluster members in the North-East and South- West seg- 
ment are likely removed by only using stars with intermediate 
brightness and masses above 4.6 M0. 

For the Padova isochrones the slope of the mass function 
decreases with increasing age going from 3 to 5 Myr from 
(-1.72 ± 0.09) to (-1.52 ± 0.09). As can be seen in Fig. HI] the 
range of initial masses along the upper part of the MS starting 
at about 20 M0 strongly decreases with age. This causes the 
same number of brighter stars in the CMD being squeezed 
into a smaller mass interval for the older ages, which therefore 
results in a flattening of the slopes with increasing age. The 
initial masses derived using the 4 Myr Padova isochrone are 
0% - 7% larger for stars with m < 37 Mq than the initial masses 
determined with the 4 Myr Geneva isochrone. For stars with 
even higher masses the Geneva isochrone yields slightly larger 
masses. This varying difference in the deduced masses causes 
the slope derived from the Geneva isochrone to be steeper by 
0.08 in comparison with the Padova isochrone of the same age. 
Nonetheless, the two derived slopes agree within the fitting 
uncertainties. 

All slopes derived binning the data into 12 bins containing 
approximately the same number of stars per bin are summa- 
rized in Table ]5] The slopes are all internally consistent: 1.) For 
each isochrone the slope derived for the South- West segment is 
flatter than the slope for the North-East segment, and the slope 
of the full sample is very close to the average of the slopes of 
both segments. 2.) Independent of the sample (NE + SW, NE, 
SW), the slope of the mass function decreases with the assumed 
cluster age for the Padova isochrones, and the use of the 4 Myr 
Geneva isochrone results in a steeper slope than for the 4 Myr 
Padova isochrone. For the North-East and South- West segment 
all slopes agree within the formal fitting uncertainties irrespec- 
tive of the isochrone. For the full sample the error margins of the 
slope derived for the 5 Myr Padova isochrone have just no over- 
lap with the error margins of the slopes derived for the 3 Myr 
Padova, the 4 Myr Padova and the 4 Myr Geneva isochrone. The 
average value of the slopes of all considered isochrones, using 
the full sample of cluster members, therefore provides a robust 
estimate for the mass function slope of the Quintuplet cluster. 
The average slope is - 1 .66, which is the same value as the slope 
of the best-fitting 4 Myr Padova isochrone. The maximum differ- 
ence between this average and the four regarded slopes is 0.14, 
which provides a conservative estimate of the uncertainty of the 
average slope. 

Our best value of the slope of the present-day mass function 
of the Quintuplet cluster for stars within a radius of 0.5 pc from 
the cluster centre, an inital mass of m,„/, > 5 M0 and excluding 
spectroscopically identified Wolf-Rayet stars, is c = -1.66 + 
0.14. It should be noted that we determined the slope using the 
initial masses as inferred from the isochrones. Furthermore, bi- 
narity is not accounted for, as we cannot observe a binary se- 
quence in the CMD of the Quintuplet cluster. Th erefore, the re- 
ported slopes refer to the system mass function. IWeidner et alJ 
(2009*) performed a numerical study to determine the influence 
of unresolved multiple systems on the initial mass function. 
Assuming 100% of the stars being part of multiple systems and 
using three different pairing methods they find that the differ- 
ence of the slopes of the single star and the observed system 
IMF for the high mass stars (m > 2 Mq) is in general smaller 
than the usual error bars of observational slopes. In general, the 
syste m IMF tends to be steeper by about 0.1 than the single star 
IMF Id a Rio et al ] (l2009f) derive a maximum difference between 
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Table 5. Overview of derived slopes of the mass function binning the data into 12 bins containing approximately the same number 
of stars. 



Isochrone name 


segment" 


No. of stars 








a 


Ao-fl," 


^^^binning 












(Mq) 


(Mq) 








3 Myr Padova 


NE + 


sw 


219 


18-19 


5.5 


65.8 


-1.72 


0.09 


0.02 


4 Myr Padova 


NE + 


sw 


257 


21-22 


4.6 


46.7 


-1.66 


0.09 


0.01 


5 Myr Padova 


NE + 


sw 


283 


24-25 


4.0 


34.5 


-1.52 


0.09 


0.02 


4 Myr Geneva 


NE + 


sw 


257 


21-22 


4.5 


47.3 


-1.74 


0.09 


0.01 


3 Myr Padova 


NE 




100 


8-9 


5.5 


65.8 


-1.75 


0.13 


0.06 


4 Myr Padova 


NE 




119 


10-11 


4.6 


46.4 


-1.68 


0.13 


0.01 


5 Myr Padova 


NE 




131 


11-12 


4.0 


34.5 


-1.60 


0.13 


0.01 


4 Myr Geneva 


NE 




119 


10-11 


4.5 


47.3 


-1.76 


0.13 


0.01 


3 Myr Padova 


SW 




119 


10-11 


5.5 


65.6 


-1.71 


0.13 


0.01 


4 Myr Padova 


sw 




138 


11-12 


4.6 


51.0 


-1.64 


0.12 


0.02 


5 Myr Padova 


SW 




152 


12-13 


4.0 


36.3 


-1.50 


0.12 


0.02 


4 Myr Geneva 


sw 




138 


11-12 


4.5 


48.1 


-1.68 


0.12 


0.02 



Notes. Segment in the proper motion diagram (Fig. [SJl. *** Number of stars per bin. Average of the slopes derived for 1000 realizations of 
randomly distributing the remainder of the division of the number of stars by the number of bins into the bins by increasing the number of stars 
in the selected bin by one. Formal uncertainty of the linear fit. Maximum difference between the slopes due to the random distribution of the 
surplus stars. 



the single star and the system IMF of 0.2 for the mass range 
of 1 < Mq < 20 by using random pairing and varying the bi- 
nary fraction between and 1. The flat MF slope of a = -1.66 
observed in the Quintuplet cluster within a radius of 0.5 pc com- 
pared to the canonical slope of -2.3 is therefore not mimicked 
by non-resolved multiple systems. 

The total mass of stars in the final cluster sample amounts 
to 4355 M0 adopting the initial masses derived from the 4 Myr 
Padova isochrone and an average mass of 58 M0 for each of the 
1 1 WR stars. Extrapolating the MF down to a minimum mass of 
0.5 M0 results in a total mass of the Quintuplet cluster within a 
radius of 0.5 pc of 5895 Mq. 

9. Discussion 

All derived slopes of the mass function in the central part of 
the Quintuplet cluster above a mass of 5 Mq are systematically 
flatter than the canonical slope of the inita l mass functio n of 
a — —2.3 + 0.7 for the same mass regime dKrout)all200lh . al- 
beit still marginally contained within its large 99% confidence 
limits. This indicates that the cluster within a radius of 0.5 pc is 
depleted of lower mass stars. 

This result is not unexpected with respect to findings in other 
Galactic young massive clusters (see Table |6]l, which expose 
signs of mass segregation by steepening slopes of their mass 
function for larger distances to the cluster centre. The mass func- 
tion of the young cluster NGC 3603 (age 1 - 2.5 Myr) exposes a 
gradually steeper slope for larger annuli from a = -1.31 within 
R < 0.15pc to -1.75 for OA < R < 0.9pc dHaravama et all 
l2008l) .UDto a maximum observed distance of 3.3 pc from the as- 
sumed cluster centre the slope remains almost constant ranging 
from -1.80 to -1.86. The global slope of -1.74 for 0.4 < m < 
20 Mq is well below the canonical IMF slope of -2.3, suggest- 
ing a top-heavy IMF for this cluster Westerlund 1, with an age of 
about 3 to 5 Myr, exposes a flattened MF with a = - 1 .6 for stars 
in the mass range of 3.4 < m < 27 Mq within R < 0.75 pc, 
which successiv ely steepens at large r radii up to a = -2.7 
for R > 2.1 pc dBrandner et al.l l2b08h . These general finding s 
were confirmed in a follow-up paper by iGennaro et alJ (1201 1[) . 
which drops the assumption of radial symmetry for the clus- 
ter and determines the mass function in a two-dimensional ap- 
proach. Their global mass function slope is with a = -2.55^'^^^^^^ 



even steeper than the canonical slope. The Arches cluster (age 
~ 2.5 Myr), is located at a projected distance to the Galactic 
centre of 26 pc, which is almost equal to the projected distance 
of 30 pc for the Quintuplet cluster (Figeret al. 1999b). Hence 
both clusters might have formed in the same environment, al- 
beit at different times, and evolved in the strong tidal field of 
the Galactic centre. The slope of the mass function of the cen- 
tral part of the Arches cluster was first determined bv lFiger et al.l 
( 199£a) and found t o be top-hea v y wit h a slope of a = - 1 .65 for 
0.1 < R < 0.35 pc. IStolte et all ( 12002 ) found a slightly steeper 
slope of a = -1.8 + 0.2 within R < 0.4 pc. Outside this radius 
the slope steepens to a - -2.70 + 0.7, indicating again mass 
segregation towards the cluster centre. The authors correct for a 
radial extinction gradient outside of 0.2 pc and use the present- 
day masses det ermined from a 2 Myr Geneva isochrone. A more 
recent study bv lEspinoza et al.l (l2009l) finds a much steeper slope 
of -2. 1 ± 0.2 for R < 0.4 pc consistent with a canonical IMF, but 
still a flattening towards the cluster co re with a = -1.88 + 0.20 
inside of R - 0.2 pc. .Espinoza et al.l (|2009) account for differ- 
ential extinction by individually dereddening the stars and in- 
fer initial masses instead of present-day masses from a 2.5 Myr 
Geneva isochrone. Espinoza et al. (2009) have also shown that 
variations in the MF slope caused by the choice of metallicity 
and a wider range of cluster ages (2.0 - 3.2 Myr) are smaller than 
the fitting uncertainties. The steeper slopes are then most likely 
a consequence of individual dereddening each star prior to the 
stellar mass estimation. This suggests that individual deredden- 
ing is one of the most crucial aspects under variable extinction 
conditions to obtain realistic MF slopes. For the Quintuplet anal- 
ysis presented above, individual dereddening was taken into ac- 
count as well, and the initial stellar masses were used to create 
the MF. In this respect our slopes of the mass function of the 
Quintuplet cluster should be directly comparable with their re- 
sults. However, as proper motions were not available, their mem- 
bership selection is solely based on a strict colour-cut leaving the 
remaining contamination by field stars unclear. 

For the Arches cluster. iKim et al.l d2006l) have quantified the 
effect of the internal cluster dynamics and the evaporation in the 
Galactic tidal field on the mass function measured within an an- 
nulus of 0. 19 - 0.35 pc. Their Fokker-Planck calculations and N- 
body simulations yield a steepening of the mass function by 0.1 
to 0.2 within the present cluster lifetime of 2.5 Myr. At an older 
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age of about 3-5 Myr, the much more dispersed appearance of 
the Quintuplet with respect to the Arches cluster suggests that 
the Quintuplet cluster is dynamically more evolved and more af- 
fected by tidal effects. While the dynamical evolution provides a 
tempting explanation for the flattened MF in the cluster centre, 
N-body simulations are required to confirm or disprove whether 
the flat PDMF of the Quintuplet cluster can be explained by dy- 
namical eflBcts alone. 



10. Summary and Outlook 

We analysed high spatial resolution H- and /Tj-band data of the 
Quintuplet cluster near to the Galactic centre obtained at the 
ESOA'LT with the NACO instrument. The cluster and the field 
star population were discerned based on the individual proper 
motions determined from two /iTj-band datasets with a time base- 
line of 5.0 yr. Remaining contaminants were removed by a sub- 
sequent colour-cut, and spectroscopically identified M,K-giants 
from the Liermann spectral catalogue (Liermann et al. 2009) 
were excluded. The slope of the present-day mass function of 
the Quintuplet cluster within R < 0.5 pc and m\aii > 5 M0 was 
derived for the first time. The inital masses of the individually 
dereddened cluster members were determined from three Padova 
MS-isochrones of solar metallicity and ages of 3, 4 and 5 Myr as 
well as a 4 Myr Geneva isochrone to study the impact of us- 
ing a different set of stellar evolution models. In order to avoid 
binning biases the mass function slopes were derived from bins 
with a uniform number of stars. This method produced robust 
results for each of the four isochrones. The derived slopes range 
from a - -1.52 to -1.74, where the Salpeter slope is -2.35. The 
slope of the best fitting isochrone (4 Myr Padova) as well as the 
mean of all slopes for the four different isochrones are found to 
be -1.66 + 0.14. 

The orbital velocity and the internal velocity dispersion of 
the cluster were recently derived and constrain the formation 
and dynamical evolution in the GC tidal field (Stolte et al., in 
prep.). Using this velocity information, N-body simulations are 
currently undertaken (Harfst et al., in prep.) to probe if the dis- 
ruptive effect of the GC tidal field can be claimed responsible for 
the shallow MF slope in the cluster core, or if a central overden- 
sity of high mass stars had to be present initially during the for- 
mation of the cluster to explain the PDMF. So far the mass func- 
tion of the Quintuplet cluster was only derived for its central part 
within R < 0.5 pc. The question if the mass function steepens in 
the outer parts of the cluster as in the Arches or Westerlund 1, or 
if the global mass function remains top-heavy as in the case of 
NGC 3603, will be addressed in a future contribution. Due to the 
high field star density affecting the outer regions of the cluster, 
measured proper motions will be even more essential to derive a 
clean cluster sample than for the cluster core. These studies will 
contribute to solving the question whether star formation pro- 
ceeds differently in the GC than in the spiral arms environment. 
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Table 6. Comparison of mass function slopes in the centres of Galactic young massive clusters. 
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Distance 


Mass Range 
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-1.66 ±0.15 
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NGC 3603 
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3.55 ±0.17 


3.4 - 27 
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-1.6 3 
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Fig. 1. VLT NACO JHR^ composite image of the Quintuplet 
cluster. Outside the dotted rectangle only H- and Ks-hwA data 
are available. The dashed circle with a radius of 500 pixel or 
0.5 pc indicates the region used for the derivation of the mass 
function (see Sect. I5.2l i. Due to bad AO correction, the 7-band 
dataset was unsuitable to perform photometry and astrometry 
and was used only for this composite image. 




magnitude magnliude 



Fig. 2. Plot of the astrometric uncertainty (left panels) and the 
photometric uncertainty (right panels) vs. the magnitude for all 
four NACO datasets. The plotted photometric uncertainty does 
only include the PSF fitting uncertainty. The dashed lines mark 
the Unearity limit of the respective dataset. 



IS function of the Quintuplet cluster 




6 8 10 12 14 16 18 20 

KSinserted [^39] 

Fig. 3. Difference of the inserted and recovered magnitudes of 
artificial stars inserted into the combined image of the A'5-band 
data in 2008 plotted in dependence of the magnitude. A high- 
order polynomial fit to the median and the standard deviation 
(multiplied by a factor of 1 .5) of the magnitude difference within 
magnitude bins of 1 mag are shown as well. If the absolute mag- 
nitude difference exceeds 0.20 mag and is larger than 1 .5 times 
the fit to the standard deviation, a recovered artificial star is re- 
jected and treated as non-recovered. The vertical dotted line indi- 
cates the maximum ^T^-band magnitude at /T, - 19 mag of stars 
to be used for the proper motion analysis. 
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8 10 12 14 16 18 20 l?"^ 46"" 16.5® 15.5® 14.5® 



Ks [mag] RA (2000) 

Fig. 4. Left panel: Recovery fractions of artificial stars inserted within the inner 500 pixel from the centre of the observed field 
plotted vs. the respective magnitude in the K,- (lower abscissa) or i/-band (upper abscissa). The position of the artificial stars were 
identical in the Kg- and //-band images and a formal colour of H - Ks - 1 .6 mag was assigned to each star. The full and dotted 
lines correspond to the recovery fractions for the /T^-band data in 2003 and 2008, respectively and the dash-dotted line shows the 
completeness in the //-band. The dashed line shows the total completeness for the stars after matching the two /Tj-band and the 
//-band datasets. Only stars with < 19 mag are used for the proper motion analysis, as indicated by the vertical dotted line. 
Right panel: ^T^-band image from the second epoch with the overplotted contours representing a completeness level of 50% for the 
labelled magnitudes. 
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Fig. 5. Left panel: Plot of the combined astrometric uncertainty from the two epochs of /iT^-band data plotted vs. the magnitude of the 
second epoch (for details of the eiTor estimation see Sect. [3]). The median and standard deviation of the astrometric uncertainty above 
the linearity limit (at 14.3 mag) of the dataset in 2003 with 20.0 s DIT were fitted by polynomials. This fit of the median (lower line) 
and the sum the median and standard deviation (upper line) are drawn in all three plots. Middle panel: Astrometric uncertainty of 
stars residing within a circle with r < 500 pixel =13 .6" around the centre of the observed field. Right panel: Astrometric uncertainty 
of stars residing outside this radius, which are excluded from further analysis. 
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Quintuplet proper motion (Ks 2003 - Ks 2008) 
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Fig. 6. Proper motion diagram of stars with Ks < 19 mag. The dashed line marks the direction parallel to the Galactic plane, the 
dotted line is oriented vertically with respect to the Galactic plane and splits the proper motion diagram into the North-East and 
the South-West segment. Stars within a radius of 2 cr as derived from the Gaussian fit in Fig. Q (right panel) around the origin are 
selected as cluster members. 
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Fig. 7. Left panel: Histogram of proper motions parallel to the Galactic plane. The strongly peaked distribution of the proper 
motions for the cluster stars around the origin overlaps with the distribution of the fields stars located predominantly in the South- 
West-segment (see also Fig|6]l. Middle panel: Histogram of proper motions vertical to the Galactic plane. Right panel: Histogram of 
the 2-dimensional proper motions located in the North-East-segment of the proper motion diagram. The histogram was fitted with a 
Gaussian function and a 2 cr cut was used as the selection criterion for cluster membership (red dashed line in all three panels). 
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Fig. 8. Left panel: Colour-magnitude diagram of cluster member candidates on the basis of their proper motions. A 4 Myr isochrone 
with solar metallicity, combined using a Padova MS-isochrone and a Pisa-FRANEC PMS-isochrone, shifted to a distance of 8 kpc 
and a foreground extinction of A^^ = 2.35 mag, is shown for reference. Stars, with fluxes exceeding the linearity limit of the 
detector are drawn as crosses in all figures throughout this paper. A small number of field stars is still present after the proper 
motion selection. These stars have similar proper motions as the cluster members, but most of them can be distinguished on the 
basis of their colours. For this purpose, a two-step colour-cut was applied for stars with H > 14 mag (see Sect.|6]for details). The 
two short-dashed vertical lines mark the first colour-cut at H - Ks = 1.3 and 2.3, rejecting blue fore- and red background stars. 
Highly reddened objects to the right of the second isochrone {Ak, = 2.88 mag, black solid line) are removed as the second step of the 
colour-cut. The sample of cluster stars after the colour-cut is plotted in red, stars rejected base d on their colour are drawn in black. 
Spectroscopically identified field giants from the spectral catalogue of iLiermann et al. I (l2009h are marked with diamonds and are 
removed from the final cluster sample. Right panel: Colour-magnitude diagram of stars classified as belonging to the field according 
to their proper motion (dots) and of stars removed from the member sample based on their colour or known spectral type (triangles). 
One star, classified as belonging to the field by its proper motion, has an O-star as (ambiguous) counterpart in the Liermann spectral 
catalogue (see Se ct. |6]l and is marked with a circle. The tilted dotted line is the line of reddening according to the extinction law by 
iNishivama et al.l (|2P09) running through the population of red clump stars from the Galactic bulge. 
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Fig. 9. Colour-magnitude diagram of the final cluster sample. Stars with counterparts in the spectral catalogue of i Liermann et alJ 
(|2009|) are flagged with symbols according to their spectral type (box: WR-stars, circle: OB-stars, stars with ambiguous identification 
are additionally marked with an X-cross). The horizontal dashed and short-dashed lines mark the initia l masses along the isoch rone 
in units of M0. The tilted dotted lines show the lines of reddening according to the extinction law bv lNishivama et aO (l2009l) and 
enframe the two regions in the CMD (shaded in grey in Fig. l9lfT0land[TTI). within which the isochrone has multiple intersections 
with the line of reddening, and consequently no unique mass can be inferred for a given star. 
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Fig. 10. Colour-magnitude diagram of all stars of the final cluster sample with proper motions residing in the North-East (left panel) 
or South- West segment (right panel) of the proper motion diagram (Fig.|6|l. See Figs. [8] and |9] for details. 
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Fig. 11. Comparison of the four isochrones used to determine stellar masses. For all shown isochrones, solar metalUcity, a distance 
to the cluster of 8 kpc and a foreground extinction of Ak^ = 2.35 mag is adopted. As in Fig. |9]and Fig. \W\ the dotted lines enframe 
regions in the CMD with ambiguous mass assignments (shaded in grey) and the initial masses are labelled along the isochrones. 
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Fig. 12. Detail of the transition region between PMS and MS in the member CMD (Fig.|9]l. The three Padova isochrones (with PMS 
parts derived from Pisa-FRANEC PMS stellar models), shifted to an extinction of Ak, = 2.35 mag to fit the MS, are shown as well. 
The stars comprising the overdensities referred to in Sect. |6]are marked with crosses. The clump of stars at the MS turn-on point is 
particularly striking, and matches excellently the predicted MS turn-on point of the 4 Myr Padova isochrone. 
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0,5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 

Fig. 13. Comparison of the mass function and derived slopes for different methods of binning the data and performing the linear fit. 
Only the fit and derived slope for the completeness corrected mass function are shown. Left panel: Mass function of the Quintuplet 
Cluster derived from the colour-magnitude diagram in Fig. |9] using a MS-Padova isochrone and a cluster age of 4Myr to transform 
magnitudes into stellar masses. Only stars with m > 4.6 M0, i.e. stars above the ambiguity region in the CMD due to the PMS, are 
used. Wolf-Rayet stars are not included in the mass function, as the large uncertainty of their mass might bias the derived slopes. 
The bin sizes are adjusted such that each of the 12 bins holds approximately the same number of stars. Right Panel: Mass function 
of the same data but distributing the stars into 12 bins of a uniform logarithmic width of 0.084 dex adopting the same lower and 
upper mass limits as in the left panel. The soUd line shows the weighted linear fit, the dotted line the unweighted fit. The slopes of 
both methods are in excellent agreement. 
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Fig. 14. Upper panels: Mass functions for the 4Myr MS-Padova isochrone using all stars (left), or only stars located in the North- 
East- (middle), or South-West-segment (right) of the proper motion diagram (Fig.|6]l Lower panels: Resulting mass function if stellar 
masses are derived from a 3 or 5 Myr Padova MS -isochrone (left and middle) or 4 Myr Geneva isochrone (right). For all shown mass 
functions the stars were distributed in 12 bins of variable width with equal numbers of stars per bin. As in Fig. [13] the Wolf-Rayet 
stars were removed and only stars above the ambiguity region in the CMD due to the PMS are included. The resulting minimum 
masses are 5.5, 4.6, 4.0, and 4.5 M©, for the 3, 4, 5 Myr Padova and the 4 Myr Geneva isochrone, respectively. 
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